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1  STATISTICAL  CONSIDERATIONS  OF 
EYE  DAMAGE 


In  considering  eye  protection  we  must  take  into  account  not  only  the 
average  energy  or  average  power  density  on  the  retina  but  also  random  vari¬ 
ation,  as  noted  by  Fried  (1981).  Random  variations  can  arise  in  at  least  two 
ways.  Most  surfaces,  either  natural  or  manmade,  are  extremely  rough  on  the 
scale  of  an  optical  wavelength.  Under  illumination  by  coherent  light,  a  wave 
reflected  from  a  typical  rough  surface  consists  of  contributions  from  many 
scattering  points  or  areas.  Images  formed  at  a  given  point  are  the  combina¬ 
tion  of  amplitudes  of  spread  functions,  each  arising  from  a  different  scattering 
point  on  the  surface.  The  spread  functions  add  with  different  phases,  giving 
complex  interference  patterns,  which  are  labeled  speckle.  The  eyes  will  form 
a  random  intensity  pattern  on  the  retina:  subjective  speckle. 

In  a  static,  nonturbulent  atmosphere,  a  speckle  pattern  propagates  through 
the  atmosphere  unchanged;  it  is  statistically  stationary.  When  turbulence  is 
present,  it  seems  to  modulate  the  brightness  of  each  speckle.  Turbulence 
gives  rise  to  a  random,  moving  patchiness  in  the  atmospheric  temperature 
field  and  thus  in  the  index  of  refraction.  The  combined  effect  of  turbulent 
eddies  is  to  form  a  screen,  which  alters  the  phases  of  the  waves  that  combined 
to  form  speckles. 

It  is  easy  to  imagine  scenarios  in  which  speckles  will  form  under  battlefield 
conditions.  For  example,  a  laser  beam  entering  a  helicopter  cockpit  will 
undergo  scattering  off  the  interior  surfaces.  Airborne  lasers  will  scatter  off 
the  ground,  buildings,  and  trees.  In  such  scenarios,  the  laser  beam  may 
undergo  several  such  scatters — the  statistical  properties  of  doubly  or  multiply 
scattered  laser  radiation  can  be  expected  to  differ  from  those  of  ordinary 
speckle. 


1.1  Statistics  of  Speckle 

Classical  speckle  theory  has  been  well-developed  in  recent  years  (Good¬ 
man,  1985),  but  much  less  work  has  been  done  on  multiply  scattered  speckle. 
If  a  surface  is  truly  rough  on  the  scale  of  a  wavelength,  the  field  associated 
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with  any  singly  linear  polarization  must  be  circular,  complex  Gaussian.  The 
intensity,  7,  follows  a  negative  exponential  distribution,  such  that  the  prob¬ 
ability  of  the  intensity  lying  between  7  and  7  +  dl  is 

Mi)  =  (lK}  i  >  o 

=  0  7  <  0, 

given  a  mean  intensity  7.  Since  the  probability  density  function  is  a  nega¬ 
tive  exponential,  the  fluctuations  about  the  mean  are  pronounced  and  the 
distribution  has  a  long  fat  tail,  as  compared  with  a  Gaussian  distribution. 
The  long  tail  implies  that  high  intensities,  several  times  the  mean,  are  not 
all  that  improbable. 

If  we  define  contrast  as  the  ratio  of  the  standard  deviation  to  the  mean, 
then  the  contrast  of  a  speckle  pattern  is  unity.  Because  of  its  high  contrast, 
speckle  is  highly  disturbing  to  observers.  This  is  particularly  true  if  fine 
detail  is  of  interest  and  speckle  results  in  a  significant  loss  in  resolution. 


The  probability,  Pj,  of  exceeding  some  threshold  intensity,  7r,  of  damage 

For  a  given  threshold  probability  Pj  and  intensity  damage  threshold  IT,  the 
mean  intensity  7  at  the  retina  must  be 


or  less,  a  result  discussed  by  Fried  (1981).  Thus,  if  Pj-  is  set  at  10-6,  7 
should  be  held  to  0.0724  lj  or  less,  in  order  to  insure  that  subjective  speckle 
is  below  threshold  at  a  probability  of  one  in  a  million.  This  probability  is 
for  a  single  speckle  scene.  In  a  battlefield  situation,  an  individual  will  be 
exposed  to  large  numbers  of  such  scenes,  as  both  the  laser  source  and  the 
individual  move  relative  to  the  rough  surface  giving  rise  to  the  speckle.  The 
probability  of  10“6  carries  into  the  far  reaches  of  the  tail  of  the  distribution. 
It  is  assumed  that  the  exponential  distribution  is  still  valid  in  this  region. 
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1.2  Statistics  of  Scattered  Speckle 


As  a  model  of  speckle  scattering,  consider  a  laser  beam  illuminating  a 
rough  surface,  which  then  scatters  onto  a  second  surface,  which  is  viewed  by 
an  observer.  Alternatively,  a  laser  beam  could  be  projected  through  a  ground- 
glass  screen  onto  a  rough  object,  which  is  then  viewed.  The  condition  for 
subjective  speckle  is  that  the  observer  is  close  enough  to  the  laser-illuminated 
object  to  resolve  clearly  the  details  of  the  objective  speckle  pattern  projected 
on  the  object  by  the  ground-glass  screen. 


The  speckle  pattern  from  the  first  scattering  has  a  nominal  instantaneous 
power  density  /'  and  is  a  random  quantity  governed  by  a  negative  exponential 
distribution  with  mean  /.  The  subjective  speckle  on  the  retina  will  also  have 
a  negative  exponential  distribution.  The  overall  power  density  on  the  retina 
will  then  have  a  probability  density  pa,  given  by 

p,(/|Z)=  fp(/|/>(/W/' 

Jo 


The  integral  can  be  evaluated  in  a  straightforward  fashion.  Instead  we 
consider  a  somewhat  more  general  problem,  in  which  the  calculation  of  the 
conditional  probability  p,  is  a  special  case.  Let  X  be  a  random  variable  with 
a  gamma  probability  density  function 


Px{x) 


1 


a(axY  1 


exp(-ax) 


a  >  0, 


with  scale  factor  a  and  index  i/;  T(i/)  is  the  gamma  function.  We  note  that 
a  negative  exponential  distribution  for  speckle  is  a  gamma  distribution  with 


q  =  /  *,  i/  =  1. 

Now  consider  a  second  random  variable  K,  with  a  gamma  density  function 
index  p  and  a  scale  factor  /?.  The  product  of  X  and  Y ,  Z, 

Z  =  XY, 
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has  a  probability  density  function 


pzW  =  /0°°wMpv(j)y 

=  r  I'"''l£lp  -  (“  +  § )  dz- 

From  the  theory  of  Bessel  functions  we  have  the  equality 


where  K\(y)  is  a  Bessel  function  of  imaginary  argument  of  the  second  kind, 
which  is  sometimes  known  as  a  Basset  or  MacDonald  function  (Watson, 
1944).  The  probability  pz{z )  is  then 

[2(«^)>]]  • 

Thus,  the  product  of  two  gamma-distributed  random  variables  gives  rise  to 
a  distribution  involving  K  functions. 

Returning  to  the  problem  at  hand,  because  we  have  a  product  of  negative 
exponentials,  the  probability  density  of  the  intensity  of  scattered  speckle  is 

1' 

p,(I\i)  =  ir'Kc  2  (i)’  . 

The  threshold  probability  for  scattered  speckle  is  then 

roo  _  _  roo  /  /  \  2 

PT  =  JhP.ii\rw=2r'fhK^2(1)  j,//. 

The  integral  K0(x)dx  is  tabulated  function  (Abramowitz  and  Stegun, 
1965).  For  a  threshold  probability  Pt  =  10~6,  the  allowable  average  intensity 
of  ordinary  speckle  is  0.0724  7j,  while  that  of  scattered  speckle  is  0.017  Ij. 
The  fatter  tail  of  the  Ka  distribution  lowers  the  allowable  average  intensity 
by  a  factor  of  4.3  for  scattered  speckle. 
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1.3  Turbulence  Modification  of  Speckle 


For  a  laser  beam  propagating  through  the  atmosphere  at  levels  of  tur¬ 
bulence  below  saturation,  the  probability  distribution  for  intensity  tends  to 
follow  either  a  lognormal  or  a  Rice-Nakagami  distribution  of  the  form 

2(0*)*' 

a 


p{x)  = 


-pexp  - 


a  +  0 


r 


where  I0(x)  is  a  Bessel  function  of  the  first  kind  with  imaginary  argument 
(Goodman,  1985).  We  assume  for  the  purpose  of  this  calculation  that  the 
Rice-Nakagami  distribution  is  the  appropriate  distribution  for  low  levels  of 
turbulence.  In  turn,  the  Rice-Nakagami  distribution  can  be  very  closely 
approximated  by  a  gamma  distribution  with  index  M.  The  index  M  of 
the  gamma  distribution  is  related  to  the  parameters  a  and  0  of  the  Rice- 
Nakagami  distribution  by 


(o  +  0) 

02  +  2a/? 


and  the  normalized  moments  of  the  gamma  distribution  are  determined  by 

in  _  r(n  +  M) 

{x)n  “  MT(M) ' 


The  probability  density  function  for  the  turbulence-modulated  spectrum 
is 

P(l\i)r  rPD(I\I')pQmdI', 

Jo 

where  pq  is  the  negative  exponential  distribution  appropriate  to  the  undis¬ 
turbed,  quiet  atmosphere,  and  p£>  is  the  probability  distribution  for  the  tur¬ 
bulent  modification  of  the  intensity  distribution.  We  thus  have  the  product 
of  two  gamma  distributions.  At  low  levels  of  turbulence,  the  M  distribu¬ 
tion  is  appropriate  (Strobehn,  Wang  and  Speck,  1975).  At  higher  levels  of 
turbulence,  the  intensity  distribution  associated  with  turbulence  alters  the 
negative  exponential  distribution  of  the  intensity  associated  with  speckle.  A 
distribution  of  the  form 


p(x) 


1  Mm 

r(M)  x 


x 


M- 1 


exp  — 


Mx 

x 
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r 


results  in  a  distribution  for  the  intensity  of 


„  /  M  \  s  .  M-l  .  _ 
p(I\I)  =  2  j  I  2  ^A/-i 


where  K\j(x)  is  the  Bessel  function  of  an  imaginary  argument  of  the  second 
kind  with  index  M .  The  moments  of  this  distribution  follow  from 

r  Ar„(x)x-'rfx  = 

Jo  y  ’  2  2 

(Watson,  1944;  p.  388).  The  normalized  moments  of  p(/|7)  are 

/"  r(n  +  M)T(1  +  M) 


(I)n 


r(M)Mn 


with  a  normalized  variance  of  1  + 


The  probability  Pj  of  exceeding  some  threshold  intensity  of  damage, 

* 


Pt = 2  It)  J,t  1  Km -■  2  (t 


dl , 


can  be  obtained  by  integrating  by  parts  and  using  the  recursion  relations 

2  M 


Km-i  —  Km+i  = 

Km- i  +  K\f+i  =  —2 


which  lead  to 


Pr  = 


2  (MITy 


M 


T(MY 


1 M 


■I<m 
x 

dI<M(x ) 
dx 


2  {MItY 

I 


Observation  of  turbulently  scattered  laser  light  leads  to  estimates  of  M 
(Holmes,  Lee  and  Kerr,  1980;  Phillips  and  Andrews,  1981).  With  no  tur¬ 
bulence,  normalized  variance  of  intensity  is  unity,  which  corresponds  to  a 
large  M .  As  turbulence  increases,  the  variance  rises  to  a  value  of  about 
1.25,  corresponding  to  M  =  8;  at  still  stronger  turbulence  the  variance  again 
decreases.  At  the  limiting  value  of  infinite  A/,  the  K  distribution  asymptot¬ 
ically  approaches  the  negative  exponential  distribution.  At  values  of  M  «8 
for  a  given  value  of  Pj,  the  allowable  mean  intensity  will  be  less  than  that 
for  ordinary  speckle,  because  of  the  long  tail  of  the  I\  distribution. 
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